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ABSTRACT 



A Nusselt type analysis was performed during laminar film condensa- 



tion on the inside of a rotating truncated cone. This analysis was 
employed to determine the condensing limit of a wickless heat pipe, 
rotating about its longitudinal axis. 

Performance characteristics including the effects of geometry, 
rotational speed, and the characteristics of fluid are given. A 
comparison is made between the condensing, boiling, sonic, and entrain 
ment limits for a given heat pipe geometry. 
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SUBSCRIPTS 
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refers 


to 


boiler 
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refers 


to 


condenser 
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to 
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m 


refers 


to 


mean value 
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to 
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I. INTRODUCTION 



A fc BACKGROUND 

In recent years considerable interest has been generated in the 
study of capillary -pumped heat pipes as a means of transferring large 
quantities of heat* At present, for most pipes using ordinary fluids, 
the limit on the heat transfer appears to be the rate that the condensate 
can be returned to the evaporator [1], The capillary-pumped heat pipes 
are primarily limited by the small pressure differentials obtainable 
across the meniscus of the liquid in the wicking, which gives rise to 
the pumping force* This differential pressure must overcome the losses 
of the liquid flowing in the wick, the losses in the vapor flow, and, 
depending on the orientation, the gravitational forces* A second limit, 
the onset of nucleate boiling in the evaporator portion of the wicking 
structure, may disrupt the capillary flow [2], [3]* Nucleate boiling will 
cause the formation of bubbles in the liquid condensate. These bubbles 
may then become trapped in the wicking structure, thus impeding the re- 
turn flow and limiting the operation of the pipe [2], [3], 

The limitations of the capillary-pumped heat pipes can be overcome 
through the use of centrifugal acceleration [4]. A pipe, rotating about 
its longitudinal axis, will generate its own gravitational field, indepen- 
dent of orientation, and thereby force the condensate to return to the 
boiler* 

A rotating heat pipe is a closed hollow shaft which contains a work- 
ing fluid. Rotation about its longitudinal axis generates a rotational 
force field. The rotational acceleration, together with a slight internal 
wall taper, pumps the working fluid from the condenser to the boiler 
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section. There is no porous wick structure as in a capillary heat 
pipe. The operation of a rotating heat pipe is shown schematically 
in Fig. 1. 

Heat Ip Heat Out 




Figure I' Cutaway View of Typical Rotating 

Heat Pipe 

As the pipe is rotated, an annulus of working fluid is formed at 
the largest internal radius of the shaft. Heat is applied to this 
portion of the pipe. As heat is applied, the fluid evaporates or 
boils, forming a vapor. This vapor is then at a higher pressure than 
that at the narrower end of the pipe (condenser). This pressure differ- 
ential forces the vapor from the boiler to the condenser section. Heat 
is removed at the condenser section, causing the vapor present at that 
section to condense on the walls of the shaft. The centrifugal force 
generated by the rotational motion of the pipe insures that the con- 
densate will return to the boiler section. Thus the rotating heat 
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pipe has a closed, two-phase heat transfer cycle, and the flow of the 
condensate is not retarded by the porous wick structure. 

The rotating heat pipe may be applied wherever there is a rotating 
part and a need to transfer heat. A few applications of the rotating 
heat pipe are the cooling of: a high speed drill, the rollers for a 
conveyer table of hot metals, the rotors of electric motors, and turbine 
blading. 

B. THESIS OBJECTIVES 

The objectives of this study were: a) to understand the operational 
principles 9 b) to develop an analytical model and, c) to predict the 
performance characteristics of a rotating heat pipe. 
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II. LIMITS OF OPERATION 



As in the capillary heat pipes, the heat transfer limits of the 
rotating heat pipe are controlled primarily by fluid dynamic considera- 
tions. [5]. The effect of the following limits will be discussed: a) 
the boiling limit, b) the entrainment limit, c) the sonic limit, and 
d) the condensing limit. 

A. BOILING LIMIT 

When a surface is heated in a stagnant pool of saturated liquid, 

three basic regimes of heat transfer are observed to occur. In the low 

range of temperature differentials , (T -T ) ,the heat flux is low and the 

w s 

principle mechanism of heat transfer is free convection. As the tempera- 
ture differential is increased, nucleate boiling becomes dominant. In 
this region, bubbles of the vapor are formed on the heating surface, 
grow in size, and eventually break off. After breaking away they rise 
to the surface and burst, releasing the vapor. The formation of the 
bubbles and their consequent rising agitates the liquid, causing convec- 
tive currents to be set up. This agitation intensifies the heat transfer 
mechanism causing an upturn in the heat flux curve and increasing the 
heat transfer coefficient significantly. As the temperature differential 
is increased further, the bubbles become so numerous, that they begin to 
coalesce, and a film of vapor can then blanket the surface. The insu- 
lating effect of this vapor blanket makes the film boiling regime un- 
desirable. A maximum heat flux exists for every boiler surface at which 
the transition occurs between nucleate boiling, and film boiling. 
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During the operation of the rotating heat pipe, the heat transfer 
coefficient, as mentioned above, will be increased due to nucleate boiling. 
A further increase in the coefficient will be experienced due to the ro- 
tational acceleration of the boiler [6], [7], [8]. The boiling limit in 
the rotating heat pipe will be reached upon the departure from nucleate 
boiling. This will occur as a sudden drop in the heat transfer coef- 
ficient which occurs when the maximum heat flux mentioned above has been 
reached. If the power input is held constant during this time, the sur- 
face temperature of the container will rapidly rise, usually causing the 
surface to melt. In the literature this phenomenon has been termed 
departure from nucleate boiling, boiling crisis, or burnout [9]. The 
boiling limit of the rotating heat pipe may be derived directly from 
the work done on conventional non-rotating boilers. 

Kutateladze has postulated that the transition from nucleate boiling 
to film boiling is purely a hydrodynamic phenomenon [10]. Based on this 
assumption he has obtained the correlation formula for predicting the 
burnout heat flux. 



where 




q 



t 



A 



b 



Pv 

h 

IS 



s 



p l 



= total heat transfer, BTU/hr 
= heat transfer area in the boiler, ft 

3 

= density of the vapor, Ibm/ft 
=■ latent heat of vaporization, BTU/lbm 

— surface tension, lbf/ft 

2 

— acceleration of gravity, ft/hr 

3 

— density of the fluid, lbm/ft 



( 1 ) 
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Kutateladze 1 s experimental data indicates the average value of K is .14, 
obtaining from a range of .13 to .19 for various surface conditions of 
horizontal wires and discs [11]. A counterflow conf igurattion of in- 
viscid vapor and liquid jets in the high nucleate boiling regime was 
proposed by Zuber, Tribus, and Westwater to predict the burnout point [12]. 
Their results for the burnout heat flux on a flat plate heater shows that 
K is equal to 77 7 24 (or. 13). Experimenters have substantiated Kutateladze 1 s 
formula and have shown that the burnout heat flux increases for rotation 
approximately as the rotational acceleration to the one fourth power [13], 
[14]. Thus the boiling limit of the rotating heat pipe may be adequately 
described by the correlation formulas of Kutateladze and Zuber [9], [15]. 
The heat flux equation (1) may be recast as a heat rate by multiplying 
through by the surface area of the rotating heat pipes boiler. This gives 
equation (2) 



B. ENTRAINMENT LIMIT 

During the operation of the rotating heat pipe, a counter current 
flow of the liquid condensate along the walls and vapor in the central 
core of the pipe is set up. A shearing off of the condensate by the 
vapor , holding up the condensate and thereby causing the vapor formation 
in the boiler to exceed the condensate return, will be defined as the 




( 2 ) 



where 



K = .13 - .19 



entrainment limit. High heat fluxes during steady state operation and 
certain critical heat fluxes during start up conditions may possibly 
generate counter current (vapor -condensate) flow rates capable of 



triggering entrainment of the condensate. In counter current flow, a 
condition closely related to the entrainment problem has been defined in 
the literature as flooding [16]. When this condition is reached, it is 
impossible to increase the flow rate of the liquid component without 
making a corresponding decrease in the flow rate of the vapor component 
and vice versa. Hewitt, Lacey, and Nicholls have shown that the flow 
contact length -is a strong parameter for the flooding transition at 
high liquid flow rates. They also have shown experimentally that as the 
absolute pressure is increased, the flow rates attainable before the 
flooding transition are higher [17]. In order to understand the state- 
ments and results of Hewitt, Lacey, and Nicholls, this discussion of 
droplet entrainment from a liquid film for co-current annular flow by 
C. B, Wallis is presented: 

#r For any given liquid flow rate the gas-liquid interface is 
characterized by the following physical appearance. At zero gas 
velocity the liquid has a smooth, mirror-like appearance as it 
moves along the duct. At a low gas velocity the surface starts 
to break up into a ripple pattern. As the gas velocity is 
further increased, the ripples start to bunch up and a wave is 
produced which appears to move at a velocity several times that 
of the entrained film. This wave, referred to as a roll wave, 
moves down the duct on top of the entrained film. At still 
higher gas velocities these roll waves grow smaller (yet more 
frequent) and droplet entrainment becomes evident. These waves 
are a precondition for the onset of entrainment and the source 
of the droplets. The wave tops are sheared off by the gas 
flowing over them. The roll waves appear to 'roll 1 in the 
direction of the gas flow." [18] 



The rotational acceleration in the rotating heat pipe will be a 
significant parameter in suppressing the formation of a rolling wave as 
described by Wallis [18]. It should act in the same manner as the in- 
creased absolute pressure delayed the onset of flooding in the experiments 
of Hewitt, Lacey, and Nicholls [17] In addition, since the rotating heat 
pipe is a sealed, self contained unit, the internal pressure may be 
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regulated to such a value as to ensure no liquid entrainment will 
occur. 

At the present time there is no satisfactory theory available to 
predict the onset of entrainment, but the experimental data of Hewitt, 
et al, can be used to predict the entrainment limit [17]. 



C. SONIC LIMIT 

As the heat flux in the rotating heat pipe is increased, it is 
possible that the flow rate of the vapor down the container may be 
limited by the attainment of a choked flow condition for the vapor. 

This choking of the vapor flow will definitely limit the amount of 
energy the vapor can transport down the pipe, thus limiting the opera- 
tion of the pipe. The heat transfer mechanism employed in the rotating 
heat pipe depends on the latent energy of the working fluid. The heat 
transfer rate in the pipe may then be expressed as 

q t = % h fg P> 



where 



rn v = mass flow rate of the vapor (lbm/hr) 

From continuity, the mass flow rate of the vapor may be expressed by 




(4) 



where 

u v = velocity of the vapor (ft/sec) , and 

2 

A = cross sectional area for flow (ft ) 

Combining equations (3) and (4) yields 

q t = ^v u v A % < 5 > 
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If the cross sectional area of the pipe is known where the choked condi- 
tion occurs,^ the sonic limitation on the heat rate may then be com- 
puted using equation (5) with the vapor velocity taken as the sonic 
velocity. 

D. CONDENSING LIMIT 

The condenser may limit the rate of return of the condensate to the 
boiler by two mechanism; a) the rate that energy may be extracted from 
the system, and b) the maximum flow rate of the condensate. At the 
present no solution exists in the literature for condensation on a 
truncated cone. The solution for condensation on a rotating cone has 
been carried out by Sparrow and Hartnett [19]. Their solution employs 
a boundary layer similarity approach. As is stated by Sparrow and 
Hartnett, the solution applies only to cones that are not too slender 
and in a region sufficiently removed from the apex of the cone. As can 
be seen in Fig. 1, the rotating heat pipe will generally employ a 
slender, truncated cone for a condenser section. If Sparrow and Hartnett's 
solution for condensation on a cone were employed, then the boundary condi- 
tion of the rotating heat pipe would not be fulfilled, namely that the 
velocity of the condensate is zero at the narrowest section of the con- 
denser. Hence, a truncated cone solution is necessary. However, since 
a similarity solution for the truncated cone would not apply in the 
region near the narrow end of the condenser, due to the boundary layer 
assumptions, the following Nusselt-type analysis was made. 



As slender pipes are anticipated (small cone angles) then as a 
good first approximation, it can be assumed that choked flow occurs at 
the exit area of the boiler. 
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m. THEORETICAL ANALYSIS OF FILM CONDENSATION 
ON THE INSIDE OF A ROTATING -TRUNCATED CONE 



In attacking the problem of condensation on the inside of a trun- 
cated cone the following assumptions are made* 

A* ASSUMPTIONS 

a ) Film condensation , not dropwise condensation, occurs in the 
condenser section* Thus the condensing liquid wets the container sur- 
face and condensation occurs at a sufficient rate to allow a continuous 
film over the surface. The assumption of filmwise condensation is felt 
to be the conservative choice as dropwise condensation has been shown 
to give higher heat transfer rates [20]. 

b) The film of condensate undergoes la minar flow . Cannon and 

Kays, using a pipe rotating about its longitudinal axis with annular 
flow, have reported, 

n Tube rotation apparently has the effect of stabilizing a 
laminar flow so that a higher through-flow Reynolds number is 
obtained before transition occurs. 11 [21]. 

c) The fluid properties are constant . Heat pipes are normally 
operated as nearly isothermal devices. This condition should be easily 
achieved in the rotating system and thus the fluid properties will very 
nearly be a constant. 

d) The subcooling of the condensate may be neglected . A small 
correction factor is normally applied to the latent heat of vaporization 
term in order to take into account the effect of subcooling. The neces- 
sity for this correction is determined by the comparison of the term 
C^AT (Cp is the specific heat at constant pressure ,BTU/lbm deg F, and 
AT is the difference between the wall temperature and the saturation 
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temperature, deg F) with the latent heat of vaporization. Normally 



C AT is much less than h 
P 



fg 



and thus may be neglected. 



e) The momentum changes through the condensate are small and 
thus there is essentially a static balance of forces. Sparrow and Gregg 
have shown that the momentum effects are negligible for ordinary fluids 
in laminar film condensation on a flat plate [22]. Assumption i) , a thin 
film, also indicates that momentum changes may be neglected. 

f) The vapor exerts no drag on the motion of the condensate; 
there is no interfacial s hear . Experiments have shown that for Prandtl 
numbers of ten and higher the effect on heat transfer of interfacial 
shear is negligible,. For Prandtl numbers near one the effect is only 
slight (within five percent) [23]. 

g) There is a linear di s tribu tion of temperature through the 
condensate film. Rohsenow has shown that for Prandtl numbers above one, 
the non-linear distribution of temperature through a thin film, due to 
energy convection, is hardly distinguishable [24] . Thus energy convec- 
tion is neglected and heat transfer occurs by conduction only. 

h) The centrifugal force is much greater than the force of 
gravity and thus gravity may be neglected. 

i) The thickness of the f i lm is much less than the radius of 
curvature of the condenser wall. 

j) The vapor spac e is essentially at one pressure , P^. 

There is no radial distribution of pressure in the vapor space and the 
axial distribution of pressure is nearly a constant [25], [26]. 

k) The density of the fluid is much greater than the density of 
the vapor . 

l) The rotating heat pipe is operating at a steady state condi- 
tion. 



m) Velocity gradients in the circumferential direction rela- 



tive to the pipe wall are negligible . 

By using the above assumptions and the coordinate system shown in 
Fig. 2, an analysis similar to Nusselt's original film condensation 
theory may be used [27]. 




Figure 2- Coordinate System ond Geometry 
for Condensing Limit 



A force balance may be taken on an infinitesimal fluid element 
in the X-direction as shown in Fig. 3. This leads to an expression 
for the momentum equation within the condensate. 
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Figure 3 Force Baloncein X-Direction on an 
Infinitesimal Fluid Element 



B. MOMENTUM EQUATION (X-direction) 

I F x = 27 rrPdy- 27 rr ( P + ff dx}d/ - r 2 - r dx + (6) 

{ r + ~ 1 /} t n r n * + Pf <J 2 it r 2 d x o y s i n f 

whe re 

P =* pressure, (lbf/ft^) 
r = radius, (ft) 

X - co-ordinate measuring distance along the surface, (ft) 

Y = co-ordinate measuring distance normal to surface, (ft) 

2 

T = shear stress, (lbf/ft ) 
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<fc = half cone angle (rad) 

2 

Pf ■ density of fluid, (lbm/ft ), and 
uj “ angular velocity, (l/hr) 

This equation reduces to 

0 = T “ — + PfUJT sm 

<3y dx r ' ^ 



( 7 ) 



C. MOMENTUM EQUATION (Y-direction) 

In like manner, using Fig. 4, a force balance in the Y-direction, 
neglecting the centrifugal head due to the vapor since p^ » p • , gives 

ZF =2-7rrPdx - 2tx' P+ ^ dy} dx - p^ r2rr cosf dxdy (8) 



This reduces to 



o - - 



ap 

d y 



p^ u> r cos 



t' 



( 9 ) 




Figure 4- Fcrce Balance in Y Direction on an 
Infinitesimal Fluid Element 
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Making use of Fig. 2, the radius at any point in the fluid 



condensate may be expressed as 

r(x,y) = R q + X sin <f> - Y cos p> 



( 10 ) 



where 

R * minimum wall radius in the condenser section (ft) 
o 

Substituting equation (10) into (9) for r(x,y) and integrating gives 

(id 



5 



8 



dy : ~Au> 2 cose p f(R+ Xsin<£ - Y cos</>) dy 
P y Y 



where 



? v = pressure of the vapor (Ibf/ft ) 
S =* film thickness (ft) 



This yields 



P= F^, + /3 f w 2 cosf/> {[r q + X sin^>] [S-y] - j cr-stp [ 8?- Y* - ] } (12) 



Now differentiating equation (12) with respect to X gives 



dP dP 
dx = dx n 



O 



r> / 

+ ^ ui cos tfr | R 0 — — + 8 sm</> +Xsin ip — 
d x 



d8 

X 



- Ysin <p - 8 cos p 



(13) 



dP 



If the result of equation (13) * ~ — # is substituted into equation 

o X 

(7) it becomes 



0 = - 



+ 8sin^> +Xsm p 
+ /> ( "sn / { R. ( + 



d8 

dx 



Y sin <£-8 cos/ £-] 



(14) 



n r t I 



dr 

Jv 



: 3 



Integrating equation (14) over the limits from Y to S and T 
then rearranging yields 



. . = - d — y { r--Y> , 2 c ' f : / r - + > 

t o , Jx 1 



- Scosc j> |?-y} + pfj sir.'/ :R 0 4 X si . ( Y 



where in equation (15) for laminar flow 



d 

T z It ■ 

' 7 (3 . 



and 



u^ = velocity of the fluid, (ft/hr) 

= viscosity of the fluid, (lbm/ft hr) 
Equation (15) may now be solved for the velocity profile, u^ , 
integration over the limits 0 to Y and 0 to u^ . 



Ux= ~ 



_LiP 

/*fd 



d8 

dx 



x <■ 2 J /', l cdx Y r d 

^ j'; ] J* y- --fp 0 + Xsincjb} {SY~ 



Y 

T 



Applying assumption j), ( vapor space at one pressure, ) 

equation (16) to be rewritten as 




.Is /, ~ r os 




r 

+ X sin 'f'l *' 'in / — cos 







>in 



cos 






d x 



Note that in equation (17) = u^x,y) since£ 



(x) 



to 0, 



(15) 

} 



by 



(16) 

} 

, allows 



(17) 
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D. CONTINUITY EQUATION 



The continuity equation of fluid dynamics requires that 

r f. f = u A 

where 

^ = mass flow rate of fluid, (Ibm/hr) 
u = average velocity of fluid, (ft/sec) 
This can be written as 



s 




(18) 



Substituting equation (17) into (18) and carrying out the integration 

u 2 

rfif = - Scos^] —■ [R Q +Xsin / - Scos*/’]} ( 19 ) 

i 

{ y [r ,+■ x - '^] - A, ' 4 r '* 4 > 

Utilizing assumption i) (thin film) allows equation (19) to become 



m, 



2 77-/1 U) Jl 
/'f 



■j [ y- Xsin / ] 




(20) 



Differentiating this result with respect to X gives 



dnrif 
d x 



2nf\_uj 



{[ 



,] .cTdS 

K 0 +Xsm/ — 




+ X sin 




( 21 ) 



E. ENERGY EQUATION 

If the film surface temperature is T and the wall is at T . 

s w 

the heat transfer by conduction to an element of surface area dA 



dq f = 



k,cos p 



{r 0 + Xsm />} In { 



FL + X sin cf> 



R 0 +X sin<£ - Sco 



} 



where 



Therefore 



dA^ = ?. rr 'j K Q t X sir /l lx 

A ^ - heat transfer area in the condenser, (ft) 

= thermal conductivity of the condensate film, 
(BTU/hr ft deg F) 

T g = saturation temperature, (deg F) 

T = inside wall temperature, (deg F) 
w 



' 1a r 



which can be written as 



,f T T '\ . 

ZTTkfCOS <± ( s w l/: 

R o hX -' j .} 

' n lR +Xsin/ -R-os / 1 
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V K C" •' /, ^ 

8 ro-.'M 

V R t X sin-n J 
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I n 



The assumption of a thin film, assumption i) , gives 






r s / 



R < X sit / 



<< I 



Therefore 



I n { I - 



R COS { 

R + X s i n / 
o 



<S r 

R r • X sin 



then 

is 

( 22 ) 



(23) 

(24) 
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( 25 ) 



Equation (24) then becomes 

dq,= -^ 2 f{^X I HvTjd* 

The rate of condensation is 



dq } = h fg <1m v 

Since the pipe is operating at a steady state condition, 
for the entire system requires that 

Ini! = in\ f 

v f 



( 26 ) 

continuity 



and applying this result to equation (26) gives 



dq,= h fa dm f 



(27) 



Equating equations (27) and (25), the mass rate of condensation per 
unit length is 



djpf _ 



S^lv'Xs^j T 



(28) 



Equating equations (28) and (21), then rearranging gives 

_Xj 

S 



- 8“ {Ho-^Xsio^] d_8 + sip ^ 



(29) 



F. FILM THICKNESS 

Now carrying out this integration (see Appendix A) for 8 
final results are 

1/4 



, the 



f 3 k f^f{ T s“T w } 


p e/ ^ , 

| ( ^ ) 


(2 Pi w hf 0 


i J 



QO) 



2 

The local heat transfer coefficient, h, (BTU/hr ft deg F) , is given 
by 




(31) 



G. MEAN HEAT TRANSFER COEFFICIENT 

The mean heat transfer coefficient over the whole surface area 
of the truncated cone is given by 



h m = 



W L$2R, 



• n + L c sm 4 } J 



L -C 

t 



2 7T 1 






A il 



i J ' 1 



(32) 



where 

L c = length along the wall of the condenser, (ft) 
Integrating equation (32) (see Appendix B) yields 



h 



m 



r 1 , I 


f2 k/".' 


{2R 0 + L c sin^} 1 


1 zT 2 j ' 4 

1 - 4 L- c 
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(33) 
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IV, DETERMINATION OF THE CONDENSING LIMIT 



The analysis of film condensation on a truncated cone may be ap- 
plied to the condenser section of the rotating heat pipe. This then 
will give an equation which will predict the condensing limit for the 
rotating heat pipe. The total heat rate in the condenser, q^, is given 
by 

1| = h m 7rL c( T s-y{ 2R o +L f" n ^ (34) 



The value of h is now substituted into equation (34) and 
m 

.3 



q f = tt 



2_ 

3 



kflf <jj~ hf j. 



T, 



/Li sm f 






1/4 8/3 8/3 i/a 

"} {[ R o +L c sin ^] " R ol 



(35) 



This then is the operating equation for predicting the limiting 
heat rate for the condenser. In the equation, however, the quantity 
T^ is not a readily controllable quantity. The inside wall temperature 
will actually adjust to a temperature corresponding to the heat rate 
(the saturation temperature) and the outside wall temperature. As the 
heat rate in the rotating heat pipe is increased, the internal pres- 
sure and saturation temperature will increase, accordingly requiring 
that the inside wall temperature adjust for this increase. Thus, the 
outside wall temperature of the condenser will be the controllable and 

measurable quantity. The inside wall temperature, T , then may be found 

w 

as a function of the outside wall temperature, T , and the heat rate 

wo 

by use of Fourier 8 a conduction law through the consenser wall, 

“V {V T J d « 



(36) 



where 



= thermal conductivity of the wall, (BTU/hr ft deg F) 
^ t = condenser wall thickness, ft 



Once again the cylindrical effects have been neglected due to the thin 
wall of thickness t. Integrating equation (36) from 0 to and 0 to L c 
yields 



A. OPERATING EQUATIONS 

The operating equation for the condensing limit is then an itera- 
tive solution between equation (35) and (37). 



At this point it should be noted that although operating equa- 
tions for the condensing limit of the pipe itself have been derived, 
as an engineering consideration, the outside cooling mechanism for 

the pipe is extremely important. Indeed, it may become the limiting 

2 

consideration. This then gives a third equation to be satisfied. 




(37) 




(35) 




(37) 




(38) 



2 

This is assuming that the film and wall thicknesses are very thin, 
and thus the areas are approximately equal. 
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where 



2 

h * coolant heat transfer coefficient (BTU/hr ft deg F) 
o 

= coolant temperature, (deg F) 
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V. PREDICTED PERFORMANCE CHARACTERISTICS 



In the actual application of the rotating heat pipe, various methods 
of cooling the condenser section are anticipated. The assumption is 
made here ^however, that the cooling mechanism employed is capable of 
setting the outside surface temperature of the condenser at a con- 
stant value, T . Equation (38) is thus eliminated from the solution, 
wo 

Rearranging the remaining condensing limit equations gives 



1/4 1/4 r 

/ 2 'T f M r 2 f jiRo+l 

vmtv I w lV t vJ; r ^ T : 



7 r 2 



8/3 

1 ^31 4 

,+L^in f J ~ Rp I 
sin <f> / 



(35) 



q t ' 2w K 





(37) 



These equations are a function of the physical properties of the working 
fluid, the thermal conductivity and the geometry of the pipe, and the 
operating conditions of rotational angular velocity and outside wall 
temperature T 

wo 

A rotating heat pipe with the physical dimensions as shown in Fig. 

5 was chosen in order to demonstrate the predicted performance of a 
given pipe. The boiler section was modeled as shown in Fig. 5 in order 
to insure that an adequate supply of working fluid would be available 
for an actual pipe. Also with this model the rotating heat pipe boiler 
is closely approximated by the rotating boiler tested at Lewis Research 
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Center [7]. The results of the Lewis tests can be used to design 
and check the operation of the rotating heat pipe boiler. 




Figure 5 Simplified 



Model of o Rofofmg 



Heal Pipe 



The outside wall temperature of the condenser, T was chosen as 

WO 

70 degrees Fahrenheit. For the saturation temperature of a parti- 
cular fluid, the fluid properties were determined. An IBM Dual 360 
digital computer was employed in iteratively solving the condensing 
limit equations (35) and (37). For a given angle, 0, and rotational 

speed, ou , the inside wall temperature, T , was varied until the 

w 

heat rates of equations (35) and (37) agreed within 10 BTU per hour. 
It should be noted that the condenser length, 10 inches, and the 
radius of the boiler exit, 1.25 inches, were held constant while the 
angle 0 was varied. 
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A. EFFECT OF FLUIDS 

The choice of the working fluid in a rotating heat pipe is dependent 
upon several physical properties of the fluid as well as the operating 
conditions desired. Various working fluids may be compared by the use 
of a dimensional parameter, N, defined as "the figure of merit." 

N ■ {M!} (39) 

The heat rate will be proportional to the figure of merit of 
various fluids for a given saturation temperature. Table I gives the 
variation in N for three fluids at a saturation temperature of 122 
degrees Fahrenheit . 

TABLE I 

Figure of Merit 



Fluid N 

Water 20 

Ethyl Alcohol 5 

Freon 12 2.4 



It should be noted that for ordinary fluids, water is the most 
desirable working fluid. This is primarily due to the combination of 
a relatively high thermal conductivity (which is cubed) coupled with a 
high latent heat of vaporization. 

The vaporization curve of a given fluid will determine the operating 
temperature range of the fluid. The heat rate attainable will then be 
determined by the choice of fluid and the operating temperature of the 
pipe as shown in Fig. 6. It should be noted that in Fig. 6 the geometry 
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of the rotating heat pipe is as shown in Fig. 5 with a constant rotation- 
al speed of 2700 RPM, and the half cone angle, 0, of 3 degrees. 

B. EFFECT OF GEOMETRY 

For the given geometric configuration shown in Fig. 5 the half cone 
angle, 0, was allowed to vary for a given operating condition and work- 
ing fluid. Figure 7 shows a plot of the heat rate versus the half cone 
angle for water at a saturation temperature of 212 degrees Fahrenheit. 

The rotational speed was allowed to vary from 300 to 3000 RPM. 

Figure 7 shows that a maximum heat rate is achieved at a particu- 
lar half cone angle for a given rotational speed (the dotted line). 

This peaking is due to the model used. The boiler exit diameter was held 
constant while the half cone angle, 0, was varied. Since L^.cos 0, the 
condenser axial length was also held constant, the minimum radius of the 

condenser, R . varied. As 0 increased R decreased and thus the heat 
9 o’ o 

rate versus the half cone angle curve has a maximum point. 

Figure 7 illustrates that for the model chosen, as the rotational 
speed of the pipe is increased, the half cone angle necessary to main- 
tain the maximum heat rate decreases. 

C. EFFECT OF ROTATIONAL SPEED 

The rotational speed versus the heat rate is shown in Fig. 8 for 
various half cone angles, 0. This curve is presented for water at 212 
degrees Fahrenheit in order to determine the most desirable operating 
speed range. One of the assumptions for the condensation on a truncated 
cone was that the centrifugal force is much greater than the gravity 
force (assumption h) . For the model shown in Fig. 5, this condition is 
estimated to occur at approximately 700 RPM. 
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D. OVERALL CAPABILITY 



A rotating heat pipe with the geometry shown in Fig. 5, operating 
at 2700 RPM, with a half cone angle of 1 degree, and using water as the 
working fluid was chosen to demonstrate the effects of the boiling, en- 
trainment, sonic, and condensing limits on the maximum heat rates at- 
tainable . 

The boiling limit curve was determined from equation (2), using 
both Kutateladze 1 s value of K equal to .19 and Zuber's value of 7T /24 
(•13), thus giving a range of heat rates for the boiling limit. Since 
the boiler surface of the rotating heat pipe may be approximated by a 
flat plate, for large radii of curvature, Zuber's predicted curve should 
be the most applicable for the model chosen. 

As mentioned earlier, no suitable theory is available to predict the 
entrainment limit, so that the experimental results of Hewitt, Lacey, 
and Nicholls have been used [17]. Their data was taken from a 1.25 inch 
diameter stationary tube in annular flow with a counter-current air- 
water contact length of 9 inches. They present flooding data for the 
vapor mass flow rate as a function of the liquid mass flow rate. The 
assumption of steady state operation in the rotating heat pipe allows 
the Hewitt, Lacey, and Nicholls data to be entered at the condition of 
equal liquid and vapor mass flow rates. The mass flow rate thus obtain- 
ed was then adjusted to agree with an average rotating heat pipe dia- 
meter of 2 inches and a counter-current steam-water flow condition. 

This adjusted mass flow rate of vapor was then used in equation (3) 
to give the entrainment limit. 

q < "\ h fg (3) 
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Equation (5) was used in computing the sonic limit. The assumption 
was made that choking occurred at the exit of the boiler section for a 
given saturation temperature in the condenser. 

The condensing limit was obtained by a solution of equations (35) 
and (37) as described on page 33. 

Figure 9 gives the plot of the rotational heat pipes heat rate as 
a function of the saturation temperature for the limits just discussed. 

For each of the limit curves shown in Fig. 9 there is an associated 
uncertainty. Even with this consideration, it is evident that the ro- 
tating heat pipe shown in Fig. 5 will be condenser limited, for a given 
saturation temperature. The upper limit on the heat rate attainable 
for a certain working fluid will generally^ - be governed by the ability 
of the container to withstand the saturation pressure. 



^It is possible that a container could withstand a certain satura- 
tion pressure* and yet not be compatible with the associated wall 
temperature . 
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VI. CONCLUSIONS 



1. The film thickness and local heat transfer coefficients for 
condensation on a rotating, truncated cone may be determined by a 
Nusselt type of analysis. The resulting film thickness and local heat 
transfer coefficient are not constant as for the rotating cone and disc, 
but vary with axial position. 

2. The predicted performance for the heat pipe geometry chosen in 
this study, and operating with water as the working fluid, is con- 
denser limited. For a different geometry, the entrainment limit may 
become the dominant limit. The sonic and boiling limits should pre- 
sent no practical limitation on the rotating heat pipe. 

3. Of the ordinary fluids, water has the best characteristics 
for use as the working fluid. 

4. The actual operation of the rotating heat pipe may be limited 
by the thermal conductance of the outside surface, (i.e., upon the 
method of cooling the rotating heat pipe). 
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VII. RECOMMENDATIONS FOR FURTHER STUDY 



A. ANALYTICAL 

1. An optimization study should be conducted to determine the 
effect of the internal geometry (boiler exit area and condenser shape 
as well as length) on the rotating heat pipe condensing limit* 

2. The restrictions of negligible interfacial shear, negligible 
momentum changes in the condensate, and constant vapor pressure should 
be removed from the present analysis. This then would require a solu- 
tion of the continuity, momentum, and energy equations for the con- 
densate as well as the vapor. The equations of the condensate and the 
vapor would be coupled through the boundary condition that the shear at 
the interface, in each phase, must be equal. 

The swirl component of the flow of condensate should also be taken 
into account. This would then require that the problem be solved in 
three-dimensions, probably through the use of spherical co-ordinates. 

3. The exact location of the choked flow condition in the rota- 
ting heat pipe should be determined in order to more accurately chart 
the sonic limit. 

4. A study should be undertaken to determine accurately the en- 
trainment limit of the rotating heat pipe for the condenser geometry 
shown in Fig. 5. 

B. EXPERIMENTAL 

1. An experimental model should be built and the condensing heat 
rate limit curve for the pipe should be experimentally verified for 
various working fluids. 
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2. A comparison of filmwise and dropwise condensation in the 
condenser section should be made, 

3. A further extension of the experimental procedure should be 
to have a condenser wall with axial grooves in it. This then would 
test the contributions of interfacial shear and swirl on the condensing 
limit. 
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APPENDIX A 



SOLUTION FOR THE FILM THICKNESS , 8 (X ) 



Equation (30) may be written in the form 



3kfM 
~T~ ZT 
P* w h 




= 3 {r c + Xsin^.} 8 +2 8 4 sm<£ - 



dx 



Let 
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f/^fL I; 



]s-]wJ 



R = R 0+ Xs,n ( /> 

u = S 4 



Therefore 



R dU + 




K 
si n 



} dR = 0 



Now let 



This gives 




K 

sin <j> 



3 dB dR 

8 B R 



Integrating this gives 




(A-l) 



(A-2) 



( A-3) 



( A-4) 
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where C is an arbitrary constant. Transforming equation (A-4) back 
into the original dimensions 



i 3 3 "sm^i = R o + Xsin<jb 



(I -s ^4 



(A-5) 



In order to determine the arbitrary constant, C, the boundary condi- 
tion § = 0 at X = 0 will be used. The constant, C, then is 



Now substituting the original value of K and the value of C into equa- 
tion (A-5), then rearranging, the final solution for S is 



The solution for condensation on a rotating cone, Sparrow and 
Hartnett’s solution, can be easily obtained by the use of the boundary 
condition that 8 = 0 at X = 0 for R q = 0. This condition would give, 
from equation (A-5), that C = 0. The solution for delta is then the 
solution obtained by Sparrow and Hartnett. 




3/8 



(A-6) 
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APPENDIX B 



SOLUTION FOR THE KEAN HEAT 
TRANSFER COEFFICIENT, h-^ 



Equation (33) may be written in the form 
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[ 2 kf pfw sm /: 


’ h fQ A 


h m* L \2h + L smc/>} 
c L o c y J 
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Working with just the integral portion of equation ( B- 1 ) 
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Equation ( B-2) then becomes 



U = R . + X s i n / 
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There fore 




( B-3) 
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It should be noted that 
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Carrying out the integration therefore yields 



Finally 



Z = 






(B-5) 



Z =• 



sm *f. 
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Substituting Z back into equation (B-l) , the average heat transfer co- 
efficient, h , is given by 



h m = 






3 2 2 

2 k fCJ h|g 
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